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We study the relation between luminosity and spin-period change in X-ray binary pulsars using long-term light 
curve obtained by the MAXI/GSC all-sky survey and pulse period data from the Fermi/GBM pulsar project. 
X-ray binaries, consisting of a highly magnetized neutron star and a stellar companion, originate X-ray emission 
according to the energy of the accretion matter onto the neutron star. The accretion matter also transfers the 
angular momentum at the Alfven radius, and then spin up the neutron star. Therefore, the X-ray luminosity 
and the spin-up rate are supposed to be well correlated. We analyzed the luminosity and period-change relation 
using the data taken by continuous monitoring of MAXI/GSC and Fermi/GBM for Be/X-ray binaries, GX 
304—1, A 0535-1-26, GRO J1008—57, KS 1947-^300, and 2S 1417—624, which occurred large outbursts in the 
last four years. We discuss the results comparing the obtained observed relation with that of the theoretical 
model by Ghosh & Lamb (1979). 


1. Introduction 

X-ray binary pulsars (XBPs) are systems consisting 
of magnetized neutron stars and mass-donating stel¬ 
lar companions. Since the neutron stars are strongly 
magnetized, the matter flows from the companion are 
dominated by the magnetic pressure inside the Alfven 
radius, and then funneled onto the magnetic poles 
along the magnetic field lines. The accretion mat¬ 
ter also transfers its angular momentum at the Alfven 
radius. Therefore, the pulsar spin-up rate and the 
mass accretion rate, i.e. the X-ray luminosity, are 
thought to be closely correlated (e.g. Ghosh & Lamb 
1979, hereafter GL79 Q)- The issue is relevant to the 
fundamental parameters of the neutron stars such as 
mass, radius, and magnetic field, as well as the XBP 
evolution scenarios. 

Be XBPs, which include Be stars extending circum- 
stellar disk around the equator, are one of the ma¬ 
jor XBP subgroups 0. They often exhibit large out¬ 
bursts lasting for about a few weeks to a few months 
mostly at around the orbital phase of the neutron-star 
periastron passage. During these outbursts, simulta¬ 
neous spin-up episodes are often observed (e.g. i). 
This is naturally explained by an increase in the accre¬ 
tion rate induced by the interaction with Be-star disk, 
and the associated transfer of the angular momentum 
to the neutron star via disc-magnetosphare coupling. 
These events give us an opportunity to study the re¬ 
lation between the luminosity and the spin-up rate 
quantitatively. 

In this paper, we present the study on the rela¬ 
tion using the long-term light curve obtained by the 
MAXI/GSC all-sky survey and the period change ob¬ 
tained from the archived results of Fermi/GBM pulsar 
project. These data, taken by the continuous moni¬ 


tor for over four years, enable us to investigate their 
time variations over the entire outburst activities in 
Be XBPs. We describe the observation in § [^ the 
analysis procedure in § [31 and then discuss about the 
obtained results in § H) 


2. Observation Data 

Since the MAXI (Monitor of All-sky X-ray Image; 
i) experiment onboard the International Space Sta¬ 
tion started in 2009 August, the GSC (Gas Slit Cam¬ 
era; i) , one of the two MAXI detectors, has been 
scanning almost the whole sky every 92-minute or¬ 
bital cycle in the 2-30 keV band. To obtain the long¬ 
term luminosity variation of Be XBPs covering the 
outbursts as well as the intermission/quiescence, we 
use archived GSC light-curve data in 2-20 keV band, 
which are processed with a standard procedure [g by 
the MAXI team and archived at MAXI web site^. 

The GBM (Gamma-ray Burst Monitor; Q) on¬ 
board the Fermi Gamma-Ray Space Telescope, is an 
all-sky instrument sensitive to X-rays and gamma-rays 
with energies between 8 keV and 40 MeV. The Fermi 
GBM pulsar project [1,0 provides results of tim¬ 
ing analysis of a number of positively detected X-ray 
pulsars, including theirpulsation periods and pulsed 
fluxes via the web site |Ill| since the in-orbit operation 
started in 2008 July. We utilized the archived pulse 
period data of Be XBPs. 

We selected hve Be XBPs, GX 304-1, A 0535-f26, 
GRO J1008-57, KS 1947-k300 and 2S 1417-624 from 
targets listed in the MAXI/GSC and the Fermi/GBM 
archive for this study, because they exhibited large 
outburst activities in the last four years and their sur¬ 
face magnetic fields are well determined by the cyclo- 
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ton resonance feature in the X-ray spectrum (execpt 
for 2S 1417—624). Table U summarized characteristic 
parameters of these Be XBPs and figure [T] shows the 
time variation of the bolometric luminosity calculated 
from MAXl/GSC 2-20 keV light curve data and that 
of the pulse period obtained from the Fermi/GBM 
pulsar data during outbursts for each source. 


3. Analysis 

Observed pulse-period variations of XBPs include 
two distinct effects, the intrinsic pulsar spin-period 
change and the orbital Doppler effect. In Be XBPs, 
both of them are supposed to correlate with the orbital 
phase. Therefore, it is not straightforward to resolve 
each component from the observed data. Although 
the pulse period data of XBPs in the Fermi/GBM 
archive are corrected for the orbital Doppler effect if 
their orbital elements are determined, the orbital ele¬ 
ments have not been known in all of the Be XBPs with 
our interests. Hence, we construct a semi-empirical 
model implementing both these effects and then fit it 
to the data, in an attempt to simultaneously deter¬ 
mine the intrinsic pulse period change and the orbital 
elements. 


3.1. Modeling of period change in XBPs 

We here employ the simple theoretical model of the 
pulsar spin-up by the mass accretion via disk, pro¬ 
posed by GL79 [l|. The model has been examined 
with X-ray data, and its validity and limits are well 
studied (e.g. [3,[i3|)- In this model, the pulsar spin-up 
rate -Pspin (s yr“^) is given by 

~ Pspin = 5.0 X 10 ^/r3g n{u}s)Si{M)P^^-^^Ljj (1) 
Si{M) = 

where /rso, PnS 6 , A 7 nS 0 , -I 45 , Pspin, A 37 are the mag¬ 
netic dipole moment of the neutron star in units of 
10^° G cm^, radius in 10® cm, mass in Mq, moment 
of inertia in 10^® g cm^, spin perid in s, luminosity in 
1037 gj.g g-i^ n{uJs) is a dimensionless torque that de¬ 
pends on the fastness parameter Wg and approximately 
constant at ~ 1.4 in slow rotating pulsars satisfying 
(PspinAa^) » 1- 

The equation [ 1 ] implies that the spin-up rate —Pgpin 
follows the luminosity L as —Pspin oc The 

power-law index 7 in a model of —Pspin oc L'^ ob¬ 
tained from the fit to the observed data sometimes 
disagreed with the theoretical value of 6/7 and favor 
the rather higher value of ^ 1.2 [1, [l3|- Besides this, 
the comparison of absolute spin-up rate with equa¬ 
tion [T] has been hampered by a large uncertainty in 


the bolometric luminosity correction, which is in turn 
due to beaming effects (e.g. [I,[l3)- We hence employ 
the spin-up model expressed by 

Pspin = ^^^37 (2) 

in which the power-law index 7 = 6/7 and a corre¬ 
lation factor, a = 1.7 X s d“^ (= ag) 

reduced from the equation [T] and typical neutron-star 
parameters of Pg = 1 , Af = 1 . 4 M 0 , 145 = 1 , are 
treated as free parameters. 

XBPs are also known to spin down during the qui¬ 
escence due to the propeller effects. The rate is much 
smaller than the spin-up during the outburst bright 
phases, but may not be negligible. We accounted its 
effect with a constant spin-down parameter, /3, added 
to i/ipin as an offset. 

By combining the spin-up and spin-down models 
above, the intrinsic pulsar-spin period Pspin (t) is ex¬ 
pressed by 

Pspin (t) = Po+ f Pspin (T)dr 
J To 

= Po+ [ 13} dr (3) 

J To 

where we set the time basis rg at the first periastron 
passage in the period under analysis and define the 
pulsation period at the time rg as Pg = Pspin (t'o) . 
The model equation [3] includes four free parameters, 
Pg, a, P, 7 and requires the luminosity data L^Tit) 
as a function of time. We calculated the luminosity 
from data of the MAXI/GSC 2-20 keV light curve in 
1 -d time bin assuming the source distance, the typical 
energy spectrum of a cutoff power law from the past 
results, and the source emission to be isotopic. 

The period modulation due to the binary orbital 
motion is calculated by using the binary elements, 
which consists of orbital period Pb, eccentricity e, pro¬ 
jected semi-major axis Ox sinf, epoch rg and argument 
uiQ of the periastron. The pulsar orbital velocity v\{t) 
along the line of sight is 


vi{t) = 


27rax sinz 

PbVT^ 


{cos {v{t) + Wg) -I- ecoswg} 


( 4 ) 


where v(t) is a parameter called ’true anomaly’ de¬ 
scribing the motion on the elliptical orbit and calcu¬ 
lated from the Kepler’s equation. The observed pulse 
period, Pohs{t), is then expressed by 

Pobs(t) Pspin(0 (^1 + . (5) 


3.2. Period-change model fit 

We applied the spin-period-change model, Pspin(^) 
in equation[31 to the Fermi/GBM archived period data 
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Table I Characteristic parameters of selected Be X-ray binary pulsars and the best-fit parameters (a, /?) used in the 
period change model. 


Target name 

-fpulse 

(s) 

-forbit 

(d) 

ax sin i 

( It-s ) 

e 

B 

( 10^2 G ) 

D 

1 ( kpc ) 

ajao 

0 

( 10"® s ) 

Ref. 

GX 304-1 

275 

132.19 

500 

0.5 

4.7 

2.0 

0.28 

2.0 

[15, 

A 0535-H26 

103 

111.10 

267 

0.47 

4.3 

2.4 

1.3 

3.6 

[17, 18] 

GRO J1008-57 

93 

249.48 

530 

0.68 

6.6 

5.8 

0.49 

2.5 

[19,^ 

KS 1947-1300 

18 

40.42 

137 

0.034 

1.1 

10 

3.2 

0.69 

[21,^ 

2S 1417-624 

17 

42.18 

188 

0.44 

* 

11 

(6.8)* 

0.0:fix 

[23, 2^ 


*: The surface magnetic field B has not been measured. It is assumed to be 2 x 10^^ G. 


for A 0535-b26, KS 1947-b300, and 2S 1417-624, in 
which the binary orbital effects were corrected with 
the known orbital elements. About GRO J1008—57, 
the orbital effects are not corrected in the archived 
data, but the orbital elements have been estimated by 
M- We thus fit the data to the period model, Pohs{t) 
in equation which includes the orbital effect, em¬ 
ploying the orbital elements given in [l^. About GX 
304—1, its orbital elements have not been measured. 
We fit the period data with the model Fobs (t) in which 
the orbital elements are floated. 

As results of many model-fit attempts, we found 
that the model is able to reproduce the data approx¬ 
imately with 7 ~ 1 in all of the five targets. We thus 
fix the parameter 7 at 6/7, predicted by GL79 Q , in 
order to concentrate on the correlation factor a, here¬ 
after. In figure [1] bottom panels, the obtained best-fit 
models with ^ = djl are superposed on the period 
data. The best-fit parameters are shown in table HI 
where the values of a are given by the ratio to that 
(= ao) predicted by GL79 [l|. 

4. Discussion 

We fitted pulse period variation of five Be XBPs 
observed with Fermi/GBM to the model implement¬ 
ing the spin-up due to the mass accretion via disk, 
expressed by Apin = based on GL79 [I| and 

the luminosity estimated from the MAXI/GSC light 
curve. The results show that the model successfully 
reproduce the data in all of the five samples. The ob¬ 
tained best-fit parameters imply that the correlation 
factor a from the luminosity to the spin-up rate 
P largely agree with ao predicted by GL79 [I|. The 
dispersion of the ratio, a/ao ^ 0.3 to 3, is naturally 
expected from the uncertainty in the bolometric lumi¬ 
nosity correction due to the beaming effect. 

However, the values of a/ao seems to have some 
tendency against the pulse period, the orbital period, 
and the eccentricity, which are suggested to have a re¬ 


lation with Be-XBP subgroups [l^ . This will become 
clearer with increasing data in the near future. 
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2010 / 01/04 GX 304-1 2013 / 11/04 2013 / 08/26 KS 1847+300 2014 / 06/22 





2009 / 09/26 GROJ1008-57 2014 / 04/03 



Figure 1: For each target of GX 304—1, A 0535-1-26, GRO J1008—57, KS 19474-300, and 2S 1417—624, time variation 
of luminosity estimated from MAXI/GSC 2-20 keV light curve data in 1-d time bin (top) and and that of pulse period 
during the outbursts obtained from Fermi/GBM pulsar data (bottom) are plotted. In the top panels, solid lines 
represent the luminosity data I/ 37 (t) used for the period-change model fit. In the bottom panel, solid and dash lines 
represent the best-fit period model and the inclusive orbital Doppler effects which have been corrected in A 0535-1-26, 
KS 1947-4300, and 2S 1417-624. 
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